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Abstract 

We have studied the one-dimensional periodic, symmetric Anderson model 
at half filling in the presence of disorder using finite-temperature quantum 
Monte Carlo techniques. We have examined for the first time the disorder 
both in hybridization between the local /-orbitals and the conduction elec- 
trons and in the local /-site energy, using a uniform distribution of width 
A. The /-orbital local magnetic moment, the uniform magnetic susceptibil- 
ity, the charge compressibility, and the nearest-neighbor magnetic correlation 
function have been calculated as a function of the disorder distribution width 
A. We find that the disorder in hybridization has a dramatic effect on the 
low-temperature magnetic properties exhibiting a non- Fermi-liquid behaviour, 
and that for the range of temperature studied the magnetic susceptibility can 
be scaled by a power law with an exponent that is in agreement with recent 
experiments. On the other hand, disorder in the local /-orbital energy level 
does not show a non-Fermi liquid behaviour. 



PACS: 71.10.Hf, 71.10.Fd, 71.27.+a, 75.20.Hr 



The last several years have witnessed dramatic growth in both the experimental [0-^ and 
theoretical p|-p!2| activity in a new class of heavy fermion materials which exhibit non- Fermi- 
liquid (NFL) behaviour in their physical properties at low temperatures. Some of the heavy 
fermion NFL metals, such as CeCue-xAu^;, have been associated with the proximity to a 
quantum critical point. [|T3],|T1| However, in several other cases of heavy-fermion compounds, 



such as Yi-^U^^Pds, ^ UCu^-xPdx, and Cei_2,Thj;RhSb, |^ NFL behaviour occurs 
only when the ^electron materials, consisting primarily of Ce or U intermetallics, have 
been alloyed with a nonmagnetic element. The NFL behaviour in these compounds is 
characterized by a linear resistivity at low T, a logarithmic low temperature divergence of 
the specific heat coefficient and a logarithmic or weak power law of the susceptibility. 

Several models have been proposed to explain these experimental results. These include 
exotic single- impurity models, such as the quadrapolar Kondo model, various multi- 
channel Kondo models 0] and the "compactified" Anderson impurity model. [0 Bhatt and 
Fischer ^ and Dobrosavljevic et al have examined the single- impurity Kondo problem 
in the presence of a random distribution of nonmagnetic impurities. However, all these 
calculations strictly apply to systems with a dilute collection of Kondo centres. The inclusion 
of lattice effects in the presence of disorder presents an additional challenge that has only 
recently started to be addressed. |To|-|l^ Tesanovic employed a slave boson technique. 



More recently, Miranda et al [|Tl|] and Chattopadhyay and Jarrell |]T2| have analyzed the 



effects of disorder on concentrated Kondo alloys in the limit of infinite spatial dimensions 
and of infinite U using the dynamical mean field theory [jl5[. However, this approach does 
not take into account the RKKY interactions between the /-sites, which are pertinent to the 
formation of the singlet ground state in the periodic Anderson model (PAM) [jl6|. Another 
mechanism proposed as a possible cause of the NFL behaviour is a disordered distribution 
of Kondo temperatures [Q]. Castro et al have attributed the NFL behaviour to the 



existence of Griffiths singularities close to a quantum critical point. These singularities arise 
from the interplay between the RKKY and Kondo interactions in the presence of magnetic 
anisotropy and disorder. 



In this Letter, we present the first investigation of the effects of disorder on the mag- 
netic properties of the one-dimensional periodic Anderson model (PAM) at half filling using 
quantum Monte Carlo techniques. The quantum Monte Carlo calculations treat both the 
RKKY and the Kondo interactions on an equal footing, and allow one to treat both the 
weak and strong disorder regime. We have studied the effects of disorder both in hybridiza- 
tion, V, between the local /-orbitals and the conduction electrons and in the local /-site 
energy, Ef, using a uniform distribution of width A. The former case corresponds to substi- 
tution in the ligand sites with the / sublattice remaining unchanged, whereas the latter case 
corresponds to doping directly on the ^sites. We present results for the /-orbital local mag- 
netic moment, the static uniform magnetic susceptibility, the charge compressibility, and 
the nearest-neighbor magnetic correlation function as a function of the disorder distribution 
width A for both types of disorder. The interplay of correlation and disorder can lead to 
different types of ground states. In the case of disorder in hybridization, we find that the 
magnetic susceptibility exhibits a NFL behaviour, and that it can be scaled at low tempera- 
tures with a power law with an exponent that is in excellent agreement with experiment. |^ 
On the other hand, disorder in the local /-site energy leads to a Fermi-liquid like behaviour. 

In the absence of disorder, the PAM is an interesting model since it exhibits various types 
of insulating states ranging from an antiferromagnetic insulator to a Kondo insulator or, if 



the f-c hybridization dominates, a simple band insulator. ||T6l It is believed that the PAM 



describes the competition between magnetic ordering and singlet formation in a number of 



the heavy fermion systems. |T6| The Hamiltonian for the one-dimensional PAM is 

H = -t E«.c.+i,. + H.c.) + E^J2< + I1 y^ift^^c,,^ + H.c.) 

+UjJ24< + Y.Elni. (1) 

Here, t is the hopping parameter in the c band, Uf is the Coulomb repulsion on the / band, 
Vi is the hybridization energy between the two bands, E( and E'^ are the energy levels of 
the local / and c band, respectively, and = fi^crfi,^^ and n^^. = cfcrCi,a are the density 
operators for the / and conduction electrons at site i with spin a. In the following we will 



set t = 1 and consider the case at half filhng {E^ = 0). 

The important thing to notice in Eq. (1) is that, unhke the usual periodic Anderson 
model, the local ^site parameters Vj or E( are taken here to be random numbers (static, 
uncorrelated disorder) distributed according to uniform distributions PiiVi) or P2{E{), re- 
spectively. Namely, 

^i(^.) = ^0(A-|^.-^o|). (2) 

P2{eI) = ^&{A-\eI -Ei\). (3) 

Here, 2A denotes the width of the uniform distribution for each type of disorder, and Vq and 
Eq are the average values of Vj and and Ef , respectively. We set E^ = -U//2 (symmetric 
case in the absence of disorder) and Vq = = 1 (intermediate parameter regime). In 
the intermediate coupling regime Monte Carlo simulations provide essentially exact results, 
whereas analytic approaches are most likely to fail or are inaccurate. 

We have employed a finite-temperature quantum Monte Carlo technique with an exact 
updating procedure. [|l8l ||19| The calculated physical observables of the finite disordered 



system depend strongly on the particular realization of the disorder. Therefore, we have 
to average all quantities over a sufficient number of disorder realizations and calculate the 
averaged expectation values, 

/ + 00 
l[dV,P,{VMA)m}), (4) 
-oo 

for the case of disorder in hybridization or 

l[dEiP2{Ei){A){{Ei}), (5) 

-oo 

for the case of disorder in the local Ej energy level. Here, (A) denotes the thermal expecta- 
tion value of the operator A for a given disorder configuration, which is calculated using the 
grand-canonical quantum Monte Carlo method developed by Fye. We have carried out 



simulations on an eight-site chain and /? = ^ < 8. We have made several checks for sixteen- 
site chains. The calculated results for the magnetic moment, the uniform / susceptibility, 



and the correlations functions agree with those of the 8-site calculations within the QMC 
statistical errors. This is due to the fact that the calculated observables are local and the 
size effects are quite small. For the case of the uniform /susceptibility the main contribution 
arises from the local on-site /-/ contribution. 

The distribution of tasks (disorder realizations) and the subsequent averaging procedure 
over disorder configurations was carried out in an IBM SP2 parallel machine using the 
Message-Passing Interface (MPI). We found that the statistical error is within 5% when 
the number of disorder configurations is above thirty. Thus, the measured quantities were 
averaged over thirty different disorder configurations. 

In the present study we present results of the effect of the two different kinds of disorder 
on the following observables: (i) The square of the /-orbital local moment cr, 

1 ^ 

^ = ]7E((«-4)'))> (6) 

i 

where N is the total number of lattice sites; (ii) The static / component of the uniform 
magnetic susceptibility x/, 



^f = ^T.[ ^^((KM -^fi(^)]HT(0) -<(0)])); (7) 

ij -'^ 

(iii) The charge compressibility k, 

and (iv) The nearest-neighbor magnetic correlation function C(i,i + 1), 

C^(^,J) = -<])), (9) 

with j =1 + 1. 

In the absence of disorder, the ground state of the one-dimensional PAM exhibits short- 
range magnetic correlations; the local /-electron spin moments are compensated by correla- 
tions with those of other /-electrons, as well as with those of the conduction electrons leading 
to a nonmagnetic ground state. [|l^] The charge gap in the PAM is due to hybridization in 



the absence of the Coulomb interaction Uf and is due to the on-site Coulomb interaction 
when Uf is larger than the band width of the conduction electrons, since the /-electron 
decouples from the conduction electrons. [^] For the parameter set used in this work the 
single-impurity Kondo temperature (T/^ ^ 0.182f/[^]^/^e~is^) is T^/t = 0.096 and the 
Kondo gap (A^ ~ Tf^e^/^pJ ) is Ax/t = 0.21. 



We first present results for the case of disorder in hybridization which couples the /- 
orbital and the conduction electrons. In Fig.l we plot the impurity magnetic moment 
a, the susceptibility Txf, the charge compressibility Tk, and the nearest-neighbor spin 
correlation function C{i,i + 1) as a function of the disorder width A at low temperatures, 
(3 = 8. We find that the average of the square of the /-orbital local magnetic moment, 
a, increases as the disorder width A increases. This result shows that the number of the 
unquenched local spins at low temperatures increases with disorder. Similarly, the low- 
temperature magnetic susceptibility Txf, which measures the "effective" magnetic moment, 
increases with increasing disorder width. The low value of the effective moment for A = 
at low temperatures indicates that the magnetic correlations act to screen the /-site moment 
leading to a singlet ground state. Turning on the disorder in hybridization results in a 
destruction of the singlet ground-state and a reduction of the spin gap between the singlet 



ground state and the lowest lying excited triplet state. [£T| Fig. 1 shows also that the 
nearest-neighbor spin correlation function C{i,i + 1) increases monotonically with A. These 
short-range / - / antiferromagnetic correlations have been found to be as important as the 



/ - c correlations in the formation of the singlet state in the PAM. ||16| , |20[| Consistent with 



the above results, we find that increase of disorder in hybridization tends to suppress these 
short-range antiferromagnetic / - / correlations, resulting in non-compensated / moments. 
On the other hand, the charge compressibility Tk, shown also in Fig. 1, is almost unchanged 
by disorder in hybridization. This suggests that the charge gap might be due to correlation 
effects rather than hybridization effects. 

In order to gain further insight into the effect of the hybridization-disorder on the mag- 
netic properties of the /-electrons, we present in Fig. 2 the / magnetic susceptibility x/ as a 



function of temperature for A = (PAM) and A = 1. In agreement with a previous study for 
the PAM (A = 0), we find that as the temperature is lowered Xf saturates to a constant 
value. On the other hand, the total uniform magnetic susceptibility which contains both the 
/ and the conduction electron components goes to zero due to the presence of the spin gap. 
This Fermi liquid behaviour in the absence of disorder signals the compensation of the 



/moment by the antiferromagnetic / - / spin correlations in addition to the already-present 
Kondo-like fluctuations from the antiferromagnetic c - / screening correlations. |lT6| , pO| On 



the other hand, in the presence of disorder, Xf exhibits a non- Fermi-liquid behaviour, di- 
vering at low temperatures. We find that Xf can be scaled at low temperatures by a power 
law, X/ ~ T~'^, with the exponent 7 = 0.66, which is comparable to recent experimental 
results. Q This result might be interpreted that the hybridization-disorder causes some of 
the local / sites to possess low Kondo temperatures and consequently to behave as nearly 
free magnetic moments, leading to a NFL behaviour. 

We next present results for the effect of disorder in the /-electron energy level, Ef, on 
the magnetic properties. In Fig. 3 we show the impurity local magnetic moment a, the 
magnetic susceptibility Txf, the charge compressibility Tk, and the nearest-neighbor / - / 
spin correlation function C{i,i + 1) as a function of the /-energy-level disorder width A at 
temperature T/t = 1/8. Note, that the range of A values in Fig. 3 is larger than that 
in Fig. 1. In the weak-disorder regime, A/t < 1, all quantities are almost unchanged by 
the disorder in Ef, in contrast to the corresponding results shown in Fig. 1 for the case 
of hybridization in V. This is due to the fact that the majority of the / sites are occupied, 
since the average /-electron energy level is below the Fermi energy {El = -U/2). On the 
other hand, in the strong-disorder regime, A/t > 1, some of the local /-energy levels are 
distributed above the Fermi energy as A increases further. This in turn causes the average 
local magnetic moment a to decrease with A because some of the ^orbitals become empty. 
The "effective moment" Txf increases slightly with A, indicating that even though a few 
diluted f spins cease to contribute to the compensation state, the singlet ground state is not 
destroyed. The nearest-neighbor spin correlation function C{i,i + 1) is suppressed because 



some of the effective /-electron orbitals are diluted by disorder. The low-temperature charge 
compressibility n increases with the disorder width in the strong-disorder regime. Strong 
disorder in the distribution of the E/s leads to a reduction of the charge gap and hence an 
increase of k. 

In conclusion, we have studied the effect of two kinds of disorder on the magnetic prop- 
erties of the one-dimensional PAM using quantum Monte Carlo simulations. We find that 
disorder in hybridization modifies considerably the low-temperature properties due to the 
presence of unquenched local moments. The low-temperature susceptibilty exhibits a NFL 
behaviour which is found to diverge with a power law with an exponent in agreement with 
experiment. ||^ On the other hand, disorder in the /site energy level results in a conventional 
Fermi liquid behaviour. 
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FIGURES 

FIG. 1. The square of the /-orbital local moment a (open squares), the /magnetic susceptibility 

Txf (solid squares), the charge susceptibility Tk (sold triangles) and the nearest-neighbor /- / spin 
correlation function C{i,i + 1) (open triangles) as a function of the hybridization-disorder width 
A at temperature T/t = 1/8 with Uf/t = l and V/t = l. 
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FIG. 2. The uniform / magnetic susceptibility as a function of temperature for the hybridiza- 
tion-disorder width A = O(triangles) and A = 1 (circles) for Uf/t = 1 and V/t = 1. For {T/t < 1), 
Xf can be scaled by the power law, X/ ~ T~'^, with 7 = 0.66. 
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FIG. 3. The square of the /-orbital local moment a (open squares), the /magnetic susceptibility 
Txf (solid squares), the charge susceptibility Tk (sold triangles) and the nearest-neighbor / - /spin 
correlation function C(i, i + 1) (open triangles) as a function of the /-energy- level disorder width 
A at temperature T/t = 1/8 with Uf/t = l and V/t = l. 
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